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When the glycosylated plant enzyme horseradish peroxidase 
(HRP) is conjugated to specific antibodies, it presents a power- 
ful tool for medical applications. The isolation and purification 
of this enzyme from plant is difficult and only gives low yields. 
However, HRP recombinantly produced in the yeast Pichia 
pastoris experiences hyperglycosylation, which impedes the 
use of this enzyme in medicine. Enzymatic and chemical 
deglycosylation are cost intensive and cumbersome and hither- 
to existing P. pastoris strain engineering approaches with the 
goal to avoid hyperglycosylation only resulted in physiological- 
ly impaired yeast strains not useful for protein production 
processes. Thus, the last resort to obtain less glycosylated re- 
combinant HRP from P. pastoris is to engineer the enzyme 
itself. In the present study, we mutated all the eight N-glycosy- 
lation sites of HRP CIA. After determination of the most suit- 
able mutation at each N-glycosylation site, we physiologically 
characterized the respective P. pastoris strains in the bioreactor 
and purified the produced HRP CIA glyco-variants. The bio- 
chemical characterization of the enzyme variants revealed 
great differences in catalytic activity and stability and allowed 
the combination of the most promising mutations to potential- 
ly give an unglycosylated, active HRP CIA variant useful for 
medical applications. Interestingly, site-directed mutagenesis 
proved to be a valuable strategy not only to reduce the overall 
glycan content of the recombinant enzyme but also to improve 
catalytic activity and stability. In the present study, we per- 
formed an integrated bioprocess covering strain generation, 
bioreactor cultivations, downstream processing and product 
characterization and present the biochemical data of the HRP 
glyco-library. 
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Introduction 

The heme-containing plant enzyme horseradish peroxidase 
(HRP; EC 1 . 1 1 . 1 .7) is a Class III peroxidase catalyzing the oxida- 
tion of various substrates (e.g., amines, aromatic phenols, indoles, 
phenolic acids and sulfonates) using hydrogen peroxide (H 2 0 2 ) 
as oxidant. Horseradish peroxidase exists in at least 19 different 
isoenzyme forms in the horseradish root (Armoracia rusticana), 
of which isoenzyme CIA is the most abundant and thus the most 
studied one (e.g., Dunford 1999; Veitch and Smith 2001; Veitch 
2004; Carlsson et al. 2005; Krainer et al. 2013; Spadiut and 
Herwig 2013). It is a 34 kDa monomeric oxidoreductase contain- 
ing a heme-group as well as two Ca 2+ -ions as prosthetic groups. 
The crystal structure of HRP CIA led to the identification of nine 
N-glycosylation sites of the Asn-X-Ser/Thr type, where X can be 
any amino acid but proline, of which eight are occupied when the 
enzyme is expressed in plant (Smith et al. 1990), which is why 
the molecular mass of HRP CIA increases from 34 to ~44 kDa 
(Veitch 2004; Spadiut and Herwig 2013). Due to glycosylation 
and the presence of both, the heme-group and disulfide bridges, 
the recombinant production and subsequent preparative purifica- 
tion of HRP has proven to be very difficult (Smith et al. 1990; 
Gajhede et al. 1997; Lavery et al. 2010; Spadiut et al. 2012), 
which is why HRP is still mainly isolated from plant (Lavery 
et al. 2010). However, HRP preparations from plant describe a 
mixture of isoenzymes, which seasonally varies in composition 
and concentration, and yields are extremely low (Jermyn 1952; 
Jermyn and Thomas 1954; Shannon et al. 1966). Since HRP is a 
versatile enzyme used in numerous, quite diverse industrial and 
medical applications, such as waste water treatment, fine chemical 
synthesis, immunoassays, biosensors and coupled enzyme assays 
(e.g., Krieg and Halbhuber 2003), the controllable recombinant 
production and subsequent efficient purification of single HRP 
isoenzymes is highly desired. Thus, we have not only investi- 
gated and improved the recombinant production of the isoenzyme 
HRP CIA with Pichia pastoris in the past few years (Dietzsch 
et al. 201 la,b; Kramer et al. 2012; Zalai et al. 2012; Spadiut et al. 
2013), but also developed an efficient downstream process for 
the hypermannosylated enzyme recombinantly derived from this 
yeast (Spadiut et al. 2012; Krainer, Pletzenauer, et al. 2013). 
However, for medical applications where HRP is conjugated to 
antibodies, like antibody-directed enzyme-prodrug therapy 
(Folkes and Wardman 2001; Wardman 2002) and medical 
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diagnostics (Romero et al. 1999; Huang 2001; Palmgren et al. 
2011; Dotsikas and Loukas 2012), the degree of glycosylation of 
HRP is of utmost importance, since not only the stability of the 
enzyme but also the conjugation with antibodies is expected to 
change with varying glycosylation. Besides, the untrimmed 
yeast-derived high-mannose containing glycosylation can trigger 
immune responses in humans (personal communication with Dr. 
Lisa Folkes; Gray Institute for Radiation Oncology and Biology, 
Department of Oncology, University of Oxford). 

Obviously, the availability of an enzyme without or at least 
with reduced surface glycosylation would solve abovemen- 
tioned problems. However, the biological role and importance 
of glycans for plant peroxidases is still not completely under- 
stood and is topic of numerous studies in glycobiology. So far, 
some studies report stabilizing effects of the glycans (Narhi 
et al. 1991; Wang et al. 1996), whereas other studies do not 
show such effects (Ehlers et al. 1992; Powell and Pain 1992). 
In 1990, Smith et al. (1990) were able to produce active and 
correctly folded HRP without any glycans in Escherichia coli. 
In a following study, Tams and Welinder analyzed the import- 
ance of the glycosylation for HRP in more detail (Tams and 
Welinder 1998). They showed that the removal of most of the 
glycans, except the Af-acetylglucosamine residues, by a mild 
chemical deglycosylation with trifluoromethanesulfonic acid 
resulted in a fully active, but less stable enzyme (Tams and 
Welinder 1995, 1998). Both studies showed that glycans on the 
surface of HRP affect the physicochemical properties of the 
enzyme but are not required for catalytic activity. However, 
chemical deglycosylation only left 60% of HRP active and also 
describes a quite cumbersome procedure (Tams and Welinder 
1995, 1998), which is why it is not a useful method to obtain 
unglycosylated HRP. Although recombinant proteins from 
P. pastoris can also be deglycosylated enzymatically (Sainz-Pastor 
et al. 2006; Kang et al. 2012), also this option has to be reconsid- 
ered, since (1) enzymatic deglycosylation is only quantitative 
when the target protein is denatured and (2) the additional endo- 
glycosidases have to be removed again to obtain pure product. 
Another way to control and reduce the complexity of native yeast- 
like glycosylation on glycoproteins secreted from P. pastoris is 
through glyco-engineering (Choi et al. 2003; Hamilton et al. 
2003; Vervecken et al. 2004; Hamilton and Gemgross 2007). A 
key event in such engineering is the knockout of the OCH1 gene, 
which initiates outer-chain elongation, leading to hypermannosy- 
lation. However, in a previous study, where we knocked out this 
gene, we observed that HRP with reduced glycan complexity pos- 
sessed hampered downstream processing and that the 
glyco-engineered P. pastoris strains generated were physiological- 
ly impaired, impeding efficient production processes (Krainer 
et al. 2013). 

Consequently, the last resort to efficiently produce HRP with a 
reduced amount of surface glycosylation is to glyco-engineer the 
enzyme itself. In a recent study, two selected Asn residues of 
HRP were mutated to Asp to analyze effects on the stability of 
the enzyme and to produce more properly folded HRP in E. coli 
(Asad, Khajeh, et al. 2011). Asad, Khajeh, et al. (2011) showed 
that introducing the mutations Asnl3Asp and Asn268Asp did 
not just affect the production of HRP in E. coli, but also 
increased the catalytic constants as well as the thermal stability. 
These results did not only underline the possibility of obtaining 



active and correctly folded HRP with reduced glycosylation but 
also showed that mutating the glycosylation sites may even have 
beneficial effects on catalytic activity and stability. 

In the present study, we generated a glyco-variant library of 
HRP CIA exchanging all the eight Asn serving as glycosylation 
sites by the structurally similar amino acids Asp, Gin or Ser. We 
did not only investigate the effects of the single mutations on 
enzyme activity and stability but also on protein purification fol- 
lowing an integrated bioprocess technology approach. After deter- 
mining the most suitable mutation at the single N-glycosylation 
sites, we physiologically characterized the respective P. pastoris 
strains in the controlled environment of a bioreactor. A two-step 
purification procedure, where both chromatography steps were 
performed in the flow-through mode, enabled us to recover puri- 
fied HRP glyco-variants for subsequent biochemical characteriza- 
tion. Based thereon, we combined the most suitable mutations to 
potentially obtain an unglycosylated, active HRP variant suitable 
for medical applications. Summarizing, we conducted an inte- 
grated bioprocess study and present the bioprocess technological 
and biochemical results for the HRP C 1 A glyco-variant library. 

Results and discussion 

Screening procedure 

Every transformation into P. pastoris CBS7435 Mut s yielded 
several dozens of transformants. We randomly picked five trans- 
formants per mutation and screened them for cell growth and pro- 
duction of active HRP in shake flasks. With only few exceptions, 
all the picked transformants produced active HRP; however, on 
average only three of five showed comparable growth and prod- 
uctivity. PCR analysis confirmed the presence of the target gene 
in the genome of P. pastoris. Although we did not analyze the 
exact number of gene integration events by real-time polymerase 
chain reaction, we assumed the integration of a comparable 
amount of gene copies into the yeast genome due to the fact that 
we always transformed the same amount of linearized vector 
DNA into the P. pastoris cells (i.e., 2 ug DNA). We had observed 
such a correlation in a previous study (Krainer et al. 2012). 
Although we did not check for the exact integration site of the 
target gene in the host genome, we ascribe the observed differ- 
ences in protein production during the screening procedure to a 
fair amount of non-homologous recombination of the transformed 
HRP CIA gene into the genome of P. pastoris (Naatsaari et al. 
2012), most likely at different loci in the chromosome, which 
consequently influences the accessibility of the transcription ma- 
chinery to the transformed gene (Krainer et al. 2012). In Table I, 
we compared the growth of the best transfonnant per mutation as 
well as the specific activity and thermal stability of the produced 
HRP CIA glyco-variants. Although diagnostic applications are 
normally not performed at 60°C, we determined the stability of 
HRP at this temperature since we observed nice differences for 
the glyco-variants at this temperature and the assay could be 
easily performed in the laboratory. Furthermore, we used the 
thermal stability as a measure for kinetic stability, as also dis- 
cussed elsewhere (Polizzi et al. 2007; Spadiut et al. 2009). We 
always included a P. pastoris strain expressing the unmutated 
wild-type (wt) HRP CIA in the screening experiments as control. 
Although we obtained slightly different results for the wt depend- 
ing on the screening round, we included the average values for 
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Table I. Results of screening experiments to identify the most suitable mutation at the single N-glycosylation sites of HRP CIA 
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growth, protein production, enzyme activity and thermal stability 
in Table 1 for comparison. Based on the determined specific activ- 
ity and thermal stability, we chose the most suitable mutation at 
the single N-glycosylation sites. Except for N13S all the pro- 
duced HRP glyco-variants showed catalytic activity; however, re- 
placement of Asn by Gin never turned out to be the most suitable 
mutation at any of the eight N-glycosylation sites. Since we did 
not measure any detectable extracellular protein content for N13S 
either, we speculate that this mutation caused problems in protein 
folding and/or secretion, a phenomenon described before (Zhu 
et al. 1998; Ito, Ishimaru, et al. 2007; Ito, Seri, et al. 2007; Zou 
et al. 2013). Interestingly, we observed significant differences in 
enzyme activity and stability depending on the introduced muta- 
tion (Table I), and identified three mutations which had been 
described before, namely at positions N13 and N268 (Asad, 
Khajeh, et al. 201 1) and N255 (Lin et al. 1999), respectively. 

Physiological strain characterization in the bioreactor 
The different P. pastoris strains carrying the respective mutated 
HRP gene were physiologically characterized in single dynamic 
batch cultivations in the controlled environment of a bioreactor. 
After exhaustion of glycerol which was indicated by an increase 
in the off-gas signal, a 0.5% (v/v) methanol adaptation pulse was 
applied which was followed by several, subsequent 1.0% (v/v) 
methanol pulses (an example for this procedure is illustrated for 
the P. pastoris strains expressing the wt enzyme and variant N57S 
in Figure 1 , while illustrations for the other strains are shown in 
Supplementary data, Figure SI). This dynamic strategy has re- 
peatedly proven to be a very efficient method to physiologically 



characterize P. pastoris strains in a fast and simple manner 
(Dietzsch et al. 201 la,b; Krainer et al. 2012; Zalai et al. 2012). 

In Table II, the determined strain characteristic parameters of 
all the strains are summarized. Apparently, the introduction of 
the respective recombinant HRP C 1 A gene had an impact on the 
physiology of the P. pastoris strains. Although the majority of 
the strains showed similar maximum specific growth rates on 
glycerol (u max g i y ) between 0.24 and 0.28 b~ , the strain carrying 
the gene HRP CIA N13D showed a nearly 1.3 -fold higher fi ms3i gly 
This effect was even more pronounced with respect to the spe- 
cific methanol uptake rate during the adaptation pulse (g s adapt)- 
Surprisingly, when we calculated the average-specific methanol 
uptake rate during the consecutive 1% (v/v) methanol pulses 
(<7s average Meoti)* we observed striking differences between all 
the strains. One can speculate that the reason for the altered 
strain physiology lies in the produced HRP glyco-variant itself, 
since it is known that N-glycosylation can influence protein 
folding and protein production and thus might affect cell physi- 
ology (Zhu et al. 1998; Ito, Ishimaru, et al. 2007; Ito, Seri, et al. 
2007; Zou et al. 2013). However, as shown in Table II, the 
amount of total extracellular protein for each strain at the end of 
the dynamic batch cultivation was basically the same, indicating 
that the single mutations did not cause significant problems in 
protein folding or secretion. We also analyzed the cell-free culti- 
vation broths on SDS gels and obtained the same pattern of 
protein bands at comparable intensity (graphs not shown). Thus, 
apparently not the mutated product but rather the locus of the re- 
spective introduced gene in the yeast genome had a significant 
influence on the methanol metabolism of the cells. This influ- 
ence is also obvious in both yield coefficients (biomass yield, 
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Fig. 1. Batch cultivation of a P. pastoris CBS7435 Mut strain carrying either the unmutated HRP CIA gene (designated as "wt") or the glyco-variant HRP CIA 
N57S. Al, batch cultivation with methanol pulses of wt; Bl, batch cultivation with methanol pulses of N57S. Solid black line, carbon dioxide evolution rate (CER); 
gray square, specific substrate uptake rate (q s ). A2, calculated yields of wt; B2, calculated yields of N57S. Gray triangle, carbon dioxide yield (Jco2/s); black dot, 
biomass yield (lx/s)- 



Table II. Strain characteristic parameters determined for recombinant P. pastoris strains harboring either the wt HRP CIA gene or a glyco-variant thereof and the 
amount of total extracellular protein at the end of cultivation 
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Yx/s; carbon dioxide yield, Yco2/s)- However, we did not investi- 
gate the exact locus of gene integration for the single strains in 
more detail. Closing C-balances for all cultivations confirm the 
accuracy of the calculated strain specific parameters. It is remark- 
able that although we used the same P. pastoris strain CBS7435 
Mut s , the same vector and basically the same gene except for 
single point mutations, we obtained physiologically diverging 
strains. This actually underlines the importance of a detailed 
physiological strain characterization using the dynamic method 
applying methanol pulses, especially if subsequent fed-batch cul- 
tivations for protein production are envisioned. 



Protein purification 

After cultivation, the respective HRP glyco-variant was purified 
from the cell-free cultivation broth using a previously reported 
two-step flow-through strategy (Spadiut et al. 2012; Krainer, 
Pletzenauer, et al. 2013). Total protein content and enzymatic 
activity were determined in the flow-through and the eluates 
and the respective recovery yield of HRP activity in percentage 
(R%) and the purification factor (PF) were calculated for each 
single purification step (Table III). 

After hydrophobic charge induction chromatography 
(HCIC), we recovered >80% of wt HRP CIA and of most 
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Table III. Results of the two-step purification approach for HRP CIA applying HCIC and AEC both operated in flow-through mode 

Variant HCIC AEC Combined Spec, activity (U mg ') 
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enzyme variants except for N158D and N198D. In agreement 
with our previous observations (Spadiut et al. 2012; Krainer 
et al. 2013), the whole activity was found in the flow-through. 
The remaining 5-20% of the enzymes did not elute from the 
column under the conditions applied, which actually proofs the 
existence of a variety of enzyme species in the cultivation broth 
differing in glycosylation and thus a varying degree of inter- 
action with the resin. Interestingly, for glyco-variant N158D 
and N198D, we only recovered 50% of the enzyme in total, 
which we also found in the flow-through. Apparently, by mutat- 
ing these two glycosylation sites and thus reducing the overall 
amount of surface N-glycosylation, the masking effect thereof 
was reduced leading to a different HCIC performance for these 
two enzyme variants. However, compared with the other glyco- 
sylation sites we could not identify a particular location of these 
two sites which could potentially explain this phenomenon 
(Supplementary data, Figure S2). By HCIC, the wt enzyme was 
purified 2-fold, whereas the success of purification varied 
between 1.3- and 3.4-fold for the different glyco-variants high- 
lighting the importance of an integrated bioprocess aspect — 
already little changes of protein properties, as the degree of 
surface N-glycosylation, might have a significant impact on fol- 
lowing unit operations. The difference in the PF for wt HRP 
CIA compared with previous studies where we achieved a PF 
of 7.0 (Krainer, Pletzenauer, et al. 2013) might be explained 
by the different cultivation strategies. In our previous study, we 
cultivated P. pastoris in shake flasks, where conditions were 
not controlled and limitations in nutrients and oxygen 
occurred, which is why cells were more sensitive to cell lysis. 
Consequently, more contaminating proteins were found in the 
cell-free cultivation broth. In the present study, we cultivated the 
different strains in a bioreactor where parameters, such as pH and 
temperature, were controlled and thus undesired cell lysis was 
reduced. Consequently, the cell-free cultivation broth contained 
less contaminating proteins. This is also obvious when looking 
at the cell-free cultivation broth before purification, which 
showed a specific activity of only 20 U mg -1 from shake flasks 
(Krainer, Pletzenauer, et al. 2013) but ~40 U mg 1 in the present 
study. 

In the subsequent anion exchange chromatography (AEC) 
step, we recovered >80% of the initially applied HRP for all 
enzyme variants (Table III). However, the amount of HRP we 
found in the flow-through vastly differed between the 
glyco-variants. For the wt enzyme and variants N57S, N158D, 
N155D and N268D, we found a comparable amount of ~75% 



of HRP in the flow-through, whereas for variants N13D, 
N186D, N198D and N214S, the recovery in the flow-through 
was only 50% or less. For N13D, we even only found 25% of 
the enzyme in the flow-through, whereas the rest was found 
in the eluate together with contaminating proteins. However, 
when we looked at the location of N13 in comparison to the other 
glycosylation sites we could not identify a particularity which 
could explain this phenomenon (Supplementary data, Figure S2). 
This again highlights the importance of the single N-glycosylation 
sites and the respective surface N-glycosylation for the phy- 
sicochemical properties of HRP and the applicability of the 
flow-through chromatography approach. With regard to AEC 
purification success, we obtained a PF of >3 for the wt enzyme 
and similar values for N57S, N214S and N268D. Although we 
only recovered 25% of the initial amount of N13D in the flow- 
through, this glyco-variant was purified nearly 8-fold. Also 
N158D, N186D and N255D were purified with great success 
(5-, 10- and more than 6-fold, respectively), whereas for variant 
N198D, a PF of only two was determined in the flow-through. 

With respect to the overall purification efficiency for the 
glyco-variant library of HRP CIA using a two-step flow-through 
approach, we observed vastly different results (Table III). The 
overall recovery of the initial amount of HRP after the two purifi- 
cation steps varied from only 20% to >70%. Also the obtained 
total PF varied immensely between 2.5 for variant N198D and 
>23 for N13D. Remarkably, these vast differences only origi- 
nated from single point mutations of HRP CIA and consequent 
changes in surface N-glycosylation. Hence, for some of these 
variants, such as for N198D, we recommend to adapt the down- 
stream strategy to obtain an enzyme preparation with higher 
purity. 

Biochemical enzyme characterization 

After protein purification, the enzyme variants were biochemically 
characterized. In Table IV, the kinetic constants for the substrates 
2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) 
and H 2 C>2 are shown. The here presented apparent K m of 1.60 
mM of the wt HRP CIA preparation for ABTS was higher than 
the previously published K m values of 0.27 and 0.18 mM for 
CIA preparations from plant and E. coli, respectively (Gilfoyle 
et al. 1996). However, in previous studies on recombinant HRP 
CIA from P. pastoris, K m values of 0.68 mM (Morawski et al. 
2000) and 1.01 (Krainer, Pletzenauer, et al. 2013) were reported. 
Apparently, yeast-derived HRP CIA preparations generally have a 
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Table IV. Kinetic constants of wt HRP CIA and the glyco-variants for the substrates ABTS and H 2 0 2 as well as thermal stability 



Variant 


ABTS 






H 2 0 2 






ty 2 (n 


K m (mM) 


K m „(Umg ') 




K m (mM) 


^maxOJmg- 1 ) 




wt 


1.60 


44.2 


27.7 


0.003 


16.3 


5433 


20.6 


N13D 


2.90 


47.2 


16.3 


0.005 


14.7 


3066 


28.9 


N57S 


2.98 


113 


38.1 


0.004 


23.7 


5378 


38.5 


N158D 


3.08 


16.3 


5.30 


0.005 


51.7 


10,342 


3.2 


N186D 


4.24 


77.4 


18.2 


0.004 


7.63 


2179 


18.8 


N198D 


1.21 


14.9 


12.3 


0.003 


19.1 


5795 


18.5 


N214S 


3.48 


41.1 


11.8 


0.004 


9.36 


2531 


6.3 


N255D 


1.72 


51.5 


29.9 


0.005 


21.6 


4506 


11.6 


N268D 


1.89 


32.5 


17.2 


0.003 


10.6 


3642 


61.9 



in) 



tendency for a lowered affinity to ABTS compared with prepara- 
tions from plant and E. coli, indicating a crucial role of glycosyla- 
tion for enzyme activity. 

As shown in Table IV, the K m values of the glyco-variants 
for ABTS were higher than for the wt, except for variant 
N198D. Mutating the N-glycosylation sites on the surface of 
enzyme HRP CIA also affected the reaction rate (Table IV). 
Summarizing, in terms of catalytic efficiency with ABTS only 
variants N57S and N255D showed slightly higher or similar 
values compared with the wt. A similar effect for glyco-variant 
N255D had already been described elsewhere (Lin et al. 1999). 
The other glyco-variants were characterized by an up to 5-fold 
reduced catalytic efficiency. Interestingly, in the work of Asad, 
Khajeh, et al. (2011) HRP variants N13D and N268D also 
showed higher catalytic activity compared with the wt. 
However, Asad et al. used a different reducing substrate as well 
as different assay conditions, which is why a direct comparison 
with the present study is not feasible. 

We observed similar trends of K m and V max for the substrate 
H2O2, as the majority of HRP glyco-variants showed a reduced 
catalytic efficiency compared with the wt enzyme (Table IV). 
In fact, N57S was the only glyco-variant showing similar or 
even higher catalytic efficiency with both substrates compared 
with the wt. The Michaelis-Menten kinetics for the wt enzyme 
and for variant N57S for both substrates is exemplarily shown 
in Figure 2, whereas illustrations for the other enzyme variants 
are shown in the Supplementary data, Figure S3. 

To potentially explain the observed effects of the mutations 
on the surface of the enzyme on the catalytic activity, we deter- 
mined the distance of the respective N-glycosylation site to the 
amino acid His 170, which is linked to the heme group in the 
active site (Figure 3). However, we were not able to identify a 
direct correlation between the distance of the N-glycosylation 
site to the active site and observed changes in catalytic activity. 
Only resolving the crystal structures of the single HRP CIA 
glyco-variants and the subsequent analysis of structure-func- 
tion relationships could potentially explain the observed effects 
of the respective mutation on the catalytic behavior. 

Finally, we also tested the enzymes for thermal stability, 
since it is known that mutating N-glycosylation sites on the 
surface of proteins might affect stability (Asad, Khajeh, et al. 
201 1). In order to investigate if the protein concentration affects 
thermal stability, as described elsewhere (Asad, Khajeh, et al. 
2011), two different concentrations, i.e., 0.01 and 0.02 mg 
mL , of the wt HRP were tested. In fact, we observed a huge 
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Fig. 2. Michaelis-Menten kinetics of the unmutated wt HRP CIA and the 
glyco-variant N57S for ABTS and H 2 0 2 . (A) Kinetics for ABTS, (B) kinetics 
for H 2 0 2 . Black dots, N57S; gray squares, wt. 

difference in the half-life time at 60°C, which we determined 
with 20.6 min for the less concentrated protein solution and 
with 121 min for the more concentrated one. Consequently, 
we normalized all the different HRP glyco-variant solutions 
to a concentration of 0.01 mgmL 1 before heat treatment. 
Interestingly, we found striking differences in the thermal sta- 
bility of the enzyme glyco-variants (Table IV). For N158D, 
N214S andN255D stability was significantly reduced, whereas 
N186D and N198D showed half-life times comparable to the 
wt. Interestingly, N13D and N57S showed a higher thermal sta- 
bility and N268D had an even 3-fold higher half-life time (x/ 2 ) 
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than the fully glycosylated enzyme. Comparable observations 
for variants N13D and N268D were also made by Asad, 
Kliajeh, et al. (2011). The differences between the determined 
half-life times in this study compared with the study of Asad 
et al., who determined the stability at 50°C, can be explained 
by differences in the assay conditions. Whereas Asad, Torabi, 
et al. (201 1) used a 200 mM phosphate buffer, we only used a 
50 mM phosphate buffer, which is known to positively affect 
HRP stability. Summarizing, it is remarkable that mutating gly- 
cosylation sites of a protein and thus reducing its overall glyco- 
sylation pattern does not only reduce protein stability but might 
also increase it. 

Glycosylation analysis. To prove the absence of surface 
glycosylation on the respective mutated N-glycosylation site, we 
exemplarily performed glycosylation analysis for the glyco-variant 
N57S by digesting the protein with either chymotrypsin or trypsin 
and subsequently analyzing the peptides by liquid-chromatography 
mass spectrometry. Furthermore, aliquots of the chymotryptic 
digest were incubated with PNGase A and released glycans were 
analyzed by mass spectrometry (Figure 4). The success of 
mutating Asn57 to a non-glycosylated Ser was confirmed by mass 
spectrometric analysis of chymotryptic peptides. We clearly see the 
absence of surface N-glycosylation on N57S (Figure 4). 

Combination of mutations 

To potentially obtain an enzyme variant without any 
N-glycosylation, we combined all the eight mutations described 
in chapter Screening procedure (hereafter called "mutant"). The 
resulting P. pastoris strain was again cultivated in a batch with 
consecutive methanol pulses for physiological strain character- 
ization (Table V). Again, we observed very different strain char- 
acteristic parameters, although we used the same P. pastoris 
strain CBS7435 Mut s , the same vector and basically the same 
gene except for eight point mutations. Closing C-balances for 



both cultivations confirm the accuracy of the calculated strain 
specific parameters. In contrast to the strain carrying the wt 
HRP CIA gene, we were not able to detect any HRP activity in 
the cell-free cultivation broth for the strain carrying the mutated 
gene. Only after ultrafiltration and 20-fold concentration of the 
cultivation broth, we were able to measure activity for the 
mutated HRP CIA glyco-variant. We concentrated the enzyme 
further and diafiltrated it before we determined the catalytic 
constants with ABTS and H2O2 as well as thermal stability 
(Table VI). As shown in Table VI, the combination of all eight 
mutations to obtain a HRP CI A variant without any 7V-glycans 
resulted in an enzyme variant with extremely reduced catalytic 
efficiency and thermal stability. Although the affinity towards 
ABTS basically remained the same, the catalytic activity was 
reduced nearly 300-fold. The effects for H 2 C>2 were even more 
severe, as K m was increased >8-fold and V mm decreased 
> 100-fold resulting in a nearly 1000-fold reduced catalytic effi- 
ciency. As judged by SDS-PAGE analysis, the size of this 
variant was significantly reduced compared with the glycosy- 
lated wt enzyme (figure not shown). However, since catalytic ac- 
tivity and stability were that low, making this variant not useful 
for medical applications, we did not analyze the surface 
N-glycosylation of this glyco-variant in more detail. 

Summarizing, enzyme engineering describes a valid approach 
to obtain active HRP CIA variants with a reduced amount of 
surface N-glycosylation. Although an enzyme variant where all 
the eight N-glycosylation sites were mutated hardly showed cata- 
lytic activity and thus does not describe a meaningful tool for 
medical applications, the here described glyco-library of HRP 
CIA describes a very useful basis for further enzyme engineer- 
ing approaches. Studies, where we only combine up to four 
mutations, namely N13D, N57S, N255D and N268D, and then 
express these variants in a P. pastoris OCH\ knockout strain 
(Krainer, Gmeiner, et al. 2013) to obtain a HRP CIA variant 
useful for targeted cancer treatment are currently ongoing. 
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Fig. 4. Verification of the mutated peptide LD 57 STTSF by MS. (A) The extracted ion chromatograms for the mass of LDSTTSF ([M+H] + 770.36 Da) in N57S HRP. 
The mutant exhibited the relevant peak, whose identity was confirmed by CID fragmentation as shown in (B). B-fragments with loss of H 2 0 are designated as b° 
fragments. 



Table V. Strain characteristic parameters determined for recombinant P. pastoris strains harboring either the wt HRP CIA gene or a variant where all the eight 
N-glycosylation sites were mutated (mutant) 



Strain /W g i y (h ') 4™ adapt [h] q s adapt (mmol g 1 h ') ^ s average MeOH l»s (C mo i C m oi) Jco2/s (C mo i C m ' ol ) C-balance Protein (mg mL ') 

(mmol g 1 h ') 


Wt 0.277 11.1 0.269 0.370 
Mutant 0.222 4.7 0.660 0.882 


0.013 0.92 0.93 0.08 
0.153 0.81 0.97 0.10 


Table VI. Kinetic constants of wt HRP CIA and the variant where all the eight N- 
well as thermal stability 


glycosylation sites were mutated (mutant) for the substrates ABTS and H 2 0 2 as 


Variant ABTS 


H 2 0 2 r./ 2 (min) 


K m (mM) F max (U mg" 1 ) V^JK m (U mg" 1 mM" 1 ) 


K m (mM) K max (U mg" 1 ) V m JK m (U mg" 1 mM" 1 ) 


Wt 1.60 44.2 27.7 
Mutant 1.44 0.15 0.10 


0.003 16.3 5433 20.6 
0.026 0.14 5.38 3.2 



Material and methods 

Chemicals 

Enzymes were purchased from Fermentas GmbH (Vienna, 
Austria). ABTS diammonium salt was obtained from Sigma- 
Aldrich Handels GmbH (Vienna, Austria). Difco™ yeast nitrogen 



base w/o amino acids (YNB), Bacto™ tryptone and Bacto™ 
yeast extract were obtained from Becton Dickinson and Company 
(Schwechat, Austria). Zeocin™ was obtained from Invitrogen 
(Vienna, Austria). D-Biotin was obtained from Fluka Chemia AG 
(St. Gallen, Switzerland). All other chemicals were purchased 
from Carl Roth GmbH & Co. KG (Karlsruhe, Germany). 
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Strain and gene 

P. pastoris CBS7435 Mut s (Dietzsch et al. 2011a,b; Krainer 
et al. 2012; Zalai et al. 2012; Krainer, Gmeiner, et al. 2013; 
Spadiut et al. 2013) and vector pPpT4_S harboring the HRP iso- 
enzyme CIA, which was codon-optimized for high-level expres- 
sion in P. pastoris (Krainer et al. 2012), were used in this study. 
The codon table described in Abad et al. (20 1 0) was applied for 
codon optimization. Secretion of HRP CIA to the cultivation 
broth was facilitated by an N-terminally fused prepro-signal se- 
quence of the S. cerevisiae alpha-factor. 

Site-directed mutagenesis 

The eight Asn, representing the glycosylation sites of HRP CIA, 
were mutated to either Asp, Gin or Ser, which are amino acids 
providing a certain structural similarity to Asn, by site-directed 
mutagenesis and subsequent digestion with Dpnl (Li and 
Wilkinson 1997). The mutagenic PCR was performed as: 98°C 
for 30 s; then 10 cycles of 98°C for 10 s, 57°C for 20 s, 72°C for 
1 min— 10 cycles of 98°C for 10 s, 60°C for 20 s, 72°C for 1 
min— 10 cycles of 98°C for 10 s, 63°C for 20 s, 72°C for 1 min; 
with a final incubation at 72°C for 10 min. Each reaction con- 
tained lx HF buffer (Fermentas), 0.1 ug of plasmid DNA, 2.5 U 
Phusion DNA polymerase (Fermentas), 1 0 uM of each dNTP 
and 5 pmol of each primer in a total volume of 50 uL. All 
primers are listed in Table VII and were purchased from 
Microsynth (Vienna, Austria). 

After PCR, the methylated template DNA was degraded by di- 
gestion with 10 U of Dpnl at 37°C for at least 3 h. The remaining 
PCR products were purified using the QIAquick PCR purification 
kit (QIAGEN, Vienna, Austria) and 5 uL of each purified PCR 
product were transformed into electro-competent E. coli TOP 10 F 
' cells. The successful introduction of the desired mutation and the 
absence of further mutations were confirmed by DNA sequencing 
(Microsynth). Transformation of ~2 (j,g Swal-linearized pPpT4_S 
plasmid DNA harboring the respective mutated HRP CIA gene 
(Supplementary data, Figure S4) into P. pastoris was done as 
described by Lin-Cereghino et al. (2005). Stable transfonnants 
were generated via homologous recombination between the linear- 
ized plasmid DNA and genomic yeast DNA. Selection of success- 
fully transformed clones was performed on yeast extract peptone 
dextrose medium (YPD; 10 g L _1 yeast extract, 20 g L~' peptone, 
20 g I/ 1 glucose, 20 g L _1 agar) supplemented with 100 mg IT 1 
Zeocin. 

Screening procedure 

Screening of five randomly picked P. pastoris transformants per 
mutation was done in 1000 mL shaking flasks. We also 
included a P. pastoris CBS7435 Mut s strain carrying the unmu- 
tated HRP CIA gene (henceforth designated as "wt") as well as 
an untransformed P. pastoris CBS7435 Mut s strain as negative 
control, resulting in a total of 1 7 shaking flasks per screening 
experiment. First the clones were cultivated in 10 mL buffered 
glycerol complex medium supplemented with 100 mgL -1 
Zeocin (BMGY_Zeo; 10 g L~ yeast extract; 20 g L _I peptone; 
3.4 g L 1 YNB w/o amino acids and ammonia sulfate, 10 g L _1 
(NH4) 2 S0 4 , 400 mgL -1 biotin; 1 gL~' glycerol; 0.1 M potas- 
sium phosphate buffer, pH 6.0) in 100 mL shaking flasks at 30°C 
and 230 rpm overnight. The next day the OD 60 o was measured, 
an appropriate aliquot of the culture was taken, and after 



Table VII. Oligonucleotide primers to mutate the eight Asn residues of the 
enzyme HRP CIA, which act as N-glycosylation sites, to either Asp, Gin or Ser 



N-site Name 



Sequence (5 '— >3') 



N13 N13D_fwd AAC TCT TGT CCT GAT GTG TCC AAC ATC 

N13Q fwd AAC TCT TGT CCT CAG GTG TCC AAC ATC 

N13S fwd AAC TCT TGT CCT^GrGTG TCC AAC ATC 

N13 rev AGG ACA AGA GTTATC GTA GAA GGT TGG AGT 

N57 N57D fwd TCC ATC TTG CTG GAC GAC ACT ACC TC 

N57Q fwd TCC ATC TTG CTG GAC CAG ACT ACC TC 

N57S_fwd TCC ATC TTG CTG GAC AGC ACT ACC TC 

N57 rev GTC CAG CAA GAT GGA AGC ATC ACA ACC 

N158 N158D_fwd C AGA AAC GTT GGT CTT GAC AGATCATCC 

N158Q_fwd C AGA AAC GTT GGT CTT CAG AGATCATCC 

N158S fwd C AGA AAC GTT GGT CTT AGC AGATCATCC 

N158 rev AAG ACC AAC GTT TCT GAA AGA GTC TTT CAA 
TTG 

N186 N186D fwd ATG GAT CGT CTG TAC GAC TTC TCT AAC AC 

Nl 86Q fwd ATG GAT CGT CTG TAC CAG TTC TCT AAC AC 

Nl 86S fwd ATG GAT CGT CTG TAC AGC TTC TCT AAC AC 

N 1 86_rev GTA CAG ACG ATC CAT GAT GAATCT ACATTG GTT 

N198 N198D_rwd CCA GAT CCT ACT CTG GAC ACC ACT TAC 

N198Q_fwd CCA GAT CCT ACT CTG CAG ACC ACT TAC 

N198S fwd CCA GAT CCT ACT CTG AGC ACC ACT TAC 

N198_rev CAG AGT AGG ATC TGG CAA ACC GG 

N214 N214D fwd CCA CTT AAC GGA GAC CTG TCT GC 

N214Q_fwd CCA CTT AAC GGA CAG CTG TCT GC 

N214S_fwd CCA CTT AAC GGA ^4 GC CTG TCT GC 

N214 rev TCC GTT AAG TGG GCA CAA ACC TC 

N255 N255D_fwd TTG TTC TCC TCT CCT GAC GCT ACT GAT 

N255Q_fwd TTG TTC TCC TCT CCT CAG GCT ACT GAT 

N255S fwd TTG TTC TCC TCT CCT^GC GCT ACT GAT 

N255 rev AGG AGA GGA GAA CAA CTC CTG GTC 

N268 N268D_fwd G AGA TCC TTC GCA GAC TCC ACT CAA 

N268Q fwd G AGA TCC TTC GCA G4G TCC ACT CAA 

N268S fwd G AGA TCC TTC GCA AGC TCC ACT CAA 

N268_rev TGC GAA GGA TCT CAC CAA TGG AAT G 

The mutation sites are depicted in italics. 

centrifugation the cells were resuspended in selective buffered 
methanol complex medium supplemented with 100 mgL -1 
Zeocin (BMMY_Zeo; 10 g L 1 yeast extract; 20 g L 1 peptone; 
3.4 g L~ YNB w/o amino acids and ammonia sulfate, 10 g L -1 
(NH 4 ) 2 S0 4 , 400 mgL -1 biotin; 0.5% methanol; 0.1 M potas- 
sium phosphate buffer, pH 6.0) to an OD 60 o of 1.0. The cells 
were again cultivated at 30°C and 230 rpm. Every day, 1% (v/v) 
methanol was pulsed to the culture and a 1 mL sample was 
taken, analyzed for optical density (OD 60 o), catalytic activity 
and protein content. The catalytic activity of HRP was mea- 
sured using an ABTS assay in a CuBiAn XC enzymatic robot 
(Bielefeld, Germany). Ten microliters of sample were mixed 
with 140 uL 1 mM ABTS solution (50 mM KH 2 P0 4 , pH 6.5). 
The reaction mixture was incubated at 37°C for 5 min before 
the reaction was started by the addition of 20 uL 0.078% H 2 0 2 
(v/v). Changes in absorbance at 415 nm were measured for 80 s 
and rates were calculated. The standard curve was prepared 
using a commercially available HRP preparation (Type VI-A; 
Sigma-Aldrich) in the range from 0.02 to 2.0 U mL~ . Protein 
concentrations were determined at 595 nm by the Bradford 
assay using the Sigma-Aldrich Protein Assay Kit with bovine 
serum albumin as standard in the range of 0.2-1.2 mg mL -1 . 

After 96 h of induction, the HRP glyco-variants in the 
cell-free supernatants were tested for thermal stability. Therefore, 
aliquots of 1 mL were incubated in a waterbath at 60°C for up to 
4h, before the samples were centrifuged (14,000 rpm; 10 min) 
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and the supernatants were analyzed for remaining HRP activity. 
These values were then compared with the initial activity before 
heat treatment. Based on activity and stability measurements, the 
most suitable mutation at a respective N-glycosylation site was 
chosen and the corresponding P. pastoris strain was physiologic- 
ally characterized in the bioreactor. Before bioreactor cultivation, 
the presence of the correctly mutated HRP gene in the respective 
P. pastoris transformant was verified by colony PCR using the 
primers AOXfwd (5'-ACTCCAACCTTCTACGATAACTC-3') 
and AOX rev (5 '-ACTGTGTCATGTGCTGACC-3 ') and subse- 
quent sequencing (Microsynth). 

Strain characterization in the bioreactor 
Culture media. Precultures were done in yeast nitrogen base 
medium with 100 mg L 1 Zeocin (YNBMZeo; 3.4 g L" 1 YNB 
w/o amino acids and ammonia sulfate, 10 g L~ (NH^SO^ 400 
mg L l biotin, 20 g I7 1 glucose, 0. 1 M potassium phosphate 
buffer, pH 6.0). Batch cultivations were performed in basal 
salt medium (26.7 mLL 1 85% phosphoric acid, 1.17 gL~' 
CaS0 4 -2H 2 0, 18.2 g L 1 K 2 S0 4 , 14.9 g L 1 MgS0 4 -7H 2 0, 4.13 
gL" 1 KOH, 40gL7' glycerol, 0.2 mL L l Antifoam Struktol 
J650, 4.35 mLT 1 PTM1, NH4OH as N-source). Trace element 
solution (PTM1) was made of 6.0 g IT 1 CuS0 4 -5H 2 0, 0.08 g IT 1 
Nal, 3.0 gL 1 MnS0 4 H 2 0, 0.2 gL 1 Na 2 Mo0 4 -2H 2 0, 0.02 g 
L" 1 H3BO3, 0.5 gL/ 1 CoCl 2 , 20.0 gL 1 ZnCl 2 , 65.0 g LT 1 
FeS0 4 -7H 2 0, 0.2 g L 1 biotin, 5 mL L 1 H 2 S0 4 . Induction was 
carried out in presence of 1 mM A-aminolevulinic acid. The 
concentration of the base NH4OH was detennined by titration 
with 0.25 M potassium hydrogen phthalate. 

Experimental procedure 

Preculture. Frozen stocks (-80°C) were cultivated in 100 mL 
YNBM in 1000 mL shake flasks at 30°C and 230 rpm. The 
grown preculture was transferred aseptically to the respective 
culture vessel. The inoculation volume was 10% of the final 
starting volume. 

Batch cultivation. Batch cultivations were carried out in either 
a 3 L or a 5 L working volume glass bioreactor (Infers, 
Bottmingen, Switzerland). Basal salt medium was sterilized in 
the bioreactor and pH was adjusted to pH 5.0 by using 
concentrated NH 4 OH solution after autoclaving. Sterile filtered 
trace elements were transferred to the reactor aseptically. 
Dissolved oxygen (d0 2 ) was measured with a sterilizable d0 2 
electrode (Visiferm™, Hamilton, Bonaduz, Switzerland). The 
pH was measured with a sterilizable electrode (Easyferm™, 
Hamilton, Bonaduz, Switzerland) and maintained constant with 
a PID controller using NFLtOH solution (1-2 M). Base con- 
sumption was determined gravimetrically. Cultivation tempera- 
ture was set to 30°C and agitation was fixed to 1200 rpm. The 
culture was aerated with 1.0 vvm dried air and off-gas of the 
culture was measured by using an infrared cell for C0 2 and a 
paramagnetic cell for 0 2 concentration (Servomax, Hyderabad, 
India). Temperature, pH, d0 2 , agitation as well as C0 2 and 0 2 in 
the off-gas were measured online and logged in a process infor- 
mation management system (Lucullus, Biospectra, Schlieren, 
Switzerland). After the complete consumption of the substrate 
glycerol, indicated by an increase of d0 2 and a drop in off-gas 
activity, the first methanol pulse of a final concentration of 0.5% 



(v/v) was conducted with methanol supplemented with 12 mL 
L 1 PTM1. Following pulses were performed with 1% metha- 
nol/PTMl (v/v) (Figure 2). For each pulse, two samples were 
taken to determine the concentrations of substrate and product, 
as well as dry cell weight to calculate specific rates and yields. 

Analysis of growth and expression parameters 
Dry cell weight was determined by centrifugation of 5 mL 
culture broth (5000 rpm, 4°C, 10 min) in a laboratory centri- 
fuge (Sigma 4K15, rotor 11156), washing the pellet twice with 
5 mL deionized water and subsequent drying at 105°C to a con- 
stant weight. 

Substrate concentrations 

Concentrations of methanol were determined in cell-free 
samples by HPLC (Vienna, Austria) equipped with an ion- 
exchange column (Supelcogel C-610H Sigma- Aldrich) and a re- 
fractive index detector (Agilent Technologies). The mobile phase 
was 0.1% H 3 P0 4 with a constant flow rate of 0.5 mL min -1 and 
the system was run isocratic at 30°C. Calibration was done by 
measuring standard points in the range from 0.1 to 10gL~ 
methanol. Measurements of biornass, product and substrate con- 
centration were executed in duplicates. 

Protein purification 

After bioreactor cultivation, the cell-free cultivation broth was 
diafiltrated for subsequent HCIC using a Centramate 500S TFF 
system (PALL, Vienna, Austria) with a 1 0 kDa MWCO mem- 
brane. The buffer was HCIC-A (500 mM NaCl, 20 mM NaOAc, 
pH 6.0) and the protein solution was concentrated to a final 
volume of 40-50 mL. All further steps of concentration and 
buffer change were performed using Amicon Ultra- 15 
Centrifugal Filter Units with 10 kDa MWCO (Merck Millipore; 
Vienna, Austria). The HCIC resin MEP HyperCel™ was 
obtained from PALL and HCIC was performed in flow-through 
mode: a column containing 25 mL of MEP HyperCel™ resin 
was equilibrated with at least 4 column volumes (CV) of buffer 
HCIC-A. Forty to 50 mL concentrated HRP solution in HCIC-A 
were loaded onto the column which was then washed with at 
least 5 CV of HCIC-A at a flow rate of 55 cm If . Then a step 
elution to 100% buffer HCIC-B (1.0 M NaCl, 20 mM NaOAc, 
pH 8.0) was performed. After elution, the column was washed 
with 5 CV 0.8 M NaOH before it was stored in EtOH 20%, 1.0 
M NaCl. During all the different steps, fractions of 1 0 mL were 
collected and analyzed for protein content and catalytic activity. 

HCIC flow-through fractions showing HRP activity were 
pooled, concentrated and rebuffered in AEC-A (50 mM Tris-HCl, 
pH 8.0) for subsequent AEC using an 8 mL CLM®-DEAE mono- 
lithic column (Krainer, Pletzenauer, et al. 2013) (BIAseparations, 
Ajdovscina, Slovenia). The column was equilibrated with 5 CVof 
AEC-A at a flow rate of 16.8 cm IT 1 . Diafiltrated post-HCIC pools 
were subsequently loaded onto the AEC column at an average 
linear flow rate of 16.8 cm IT 1 before a post-load wash with 5 CV 
of AEC-A was performed. Elution was performed in a single step 
from 0 to 100% AEC-B (50 mM Tris-HCl, 1.0 M NaCl, pH 8.0). 
The column was washed with 5 CVof a 1 M NaOH/1 M NaCl so- 
lution at an average linear flow rate of 33.6 cm IT 1 for column re- 
covery, before the column was stored in 20% EtOH. 
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The efficiency of each purification step was evaluated by de- 
termining the PF and the recovery yield of HRP activity in per- 
centage (R%). PF and R% were calculated by Eqs (1) and (2) 
(Krainer, Pletzenauer, et al. 2013). The suffixes "pre" and "post" 
indicate the respective values before and after a purification step. 



PF = 



specific activity 
specific activity 



post 



(1) 



volumetric activity n „ t x volume mst 

B% = 100 x ■ P — . (2) 

volumetric activity pre x volume 

pre 

Finally, the pooled active fractions after AEC were diafiltrated 
in 50 mM potassium phosphate buffer, pH 6.5, and concentrated 
to a volume of ~ 1 .5 mL for the subsequent biochemical 
characterization. 



Biochemical enzyme characterization 

Biochemical characterization of the purified HRP glyco- variants 
included the determination of the basic kinetic parameters K m 
and F max for the two substrates H2O2 and ABTS in a spectropho- 
tometer UV-1601 from Shimadzu (Korneuburg, Austria). The re- 
action mixture with a final volume of 1 .0 mL contained 20 uL of 
HRP glyco-variant, 50 mM potassium phosphate buffer, pH 6.5, 
and either varying concentrations of ABTS (0.01-10 mM) and a 
saturating concentration of H2O2 of 1.0 mM or varying concen- 
trations of H 2 0 2 (0.001-1.0 mM) and a saturating concentration 
of ABTS of 10.0 mM, respectively. The increase in absorption 
was followed at 420 nm at 30°C for 180 s. Absorption curves 
were recorded with a software program (UVPC Optional 
Kinetics software, Shimadzu). The maximum reaction rate 
(F max ) and the Michaelis constant (K m ) were calculated with the 
Sigma Plot software (Version 1 1 .0; Systat Software Inc.). 

The thermal stability of individual HRP glyco-variants was 
tested at 60°C. The residual activity towards ABTS was mea- 
sured after 1, 5, 10, 15, 30, 45, 60, 90 and 120 min of incubation 
at 60°C in a water bath. Protein concentrations were normalized 
to 0.01 mg mL to limit possible effects of the different protein 
concentrations on thermal stability (Asad, Khajeh, et al. 2011) 
and to obtain comparability. Residual activities were plotted 
versus the incubation time and the half-life times of thermal in- 
activation at 60°C (ty 2 ) were calculated using Eq. (3): 



T l/2 



In 2 
~k~' 



(3) 



&i n rate of inactivation (slope of the logarithmic residual 
activity). 

Glycosylation analysis. Purified HRP sample was buffered in 
0.1 M NH4HCO3 and reduced with dithiothreitol (5 mM) for 
45 min at 56°C and alkylated using iodoacetamide (25 mM) at 
room temperature for 30 min. The protein was precipitated with 
4 volumes of acetone for 45 min at -20°C, dried in a vacuum 
centrifuge and resuspended in 0. 1 M NH4HCO3 buffer to yield 
a protein concentration of ~1 (ig uL -1 . Digests were performed 
overnight with either chymotrypsin or trypsin (Promega, 
Mannheim, Germany) at 37°C at an enzyme-to-substrate 



ratio of 1:50 (w/w). The digested peptides were analyzed by 
liquid-chromatography mass spectrometry as follows: 1 \ig of 
sample was loaded on a BioBasic-18 column (150 x 0.32 mm; 
5 urn; Thermo Scientific, Vienna, Austria) and eluted with a 
gradient from 1 to 60% acetonitrile in 0.3% formic acid 
buffered to pH 3.0 at flow rate of 6 uLmin. Eluted peptides 
were analyzed on an Ultima Global Q-TOF mass spectrometer 
(Waters, Manchester, UK) operated in positive-ion mode, which 
was previously calibrated with a cesium iodide standard in the 
range of 400-1800 m/z. Additionally, the peptide harboring the 
site N57S within the mutated HRP variant was subjected to 
collision-induced dissociation MS-MS with Argon as collision 
gas. Data were manually evaluated and deconvoluted using the 
Software MassLynx V4.00.00 (Waters). 

An aliquot of the chymotryptic digest was heat inactivated 
and then incubated with 0.03 mU PNGase A (Proglycan, 
Vienna, Austria) in 50 mM citrate buffer, pH 5.5. Glycans were 
purified using porous graphitic carbon cartridges (Thermo 
Scientific) as described (Pabst et al. 2012). Glycans were ana- 
lyzed by mass spectrometry as described in chapter 
Glycosylation analysis for peptides with the sole divergence of 
using a 100 x 0.32 mm hypercarb column (Thermo Scientific) 
and 1 h gradient from 1 to 50% acetonitrile. 

Supplementary Data 

Supplementary data for this article are available online at http:// 
glycob.oxfordj ournals .org/. 
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ABTS, 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid); 
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carbon dioxide yield (C-molC-mol 1 ). 
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